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A convincing identification of dark matter (DM) particles can probably be achieved only through 
a combined analysis of different detections strategies, which provides an effective way of removing 
degeneracies in the parameter space of DM models. In practice, however, this program is made 
complicated by the fact that different strategies depend on different physical quantities, or on the 
same quantities but in a different way, making the treatment of systematic errors rather tricky. We 
discuss here the uncertainties on the recoil rate in direct detection experiments and on the muon 
rate induced by neutrinos from dark matter annihilations in the Sun, and we show that, contrarily 
to the local DM density or overall cross section scale, irreducible astrophysical uncertainties affect 
the two rates in a different fashion, therefore limiting our ability to reconstruct the parameters of the 
dark matter particle. By varying within their respective errors astrophysical parameters such as the 
escape velocity and the velocity dispersion of dark matter particles, we show that the uncertainty 
on the relative strength of the neutrino and direct-detection signal is as large as a factor of two for 
typical values of the parameters, but can be even larger in some circumstances. 



PACS numbers: 95.35.+d, 14.60.Lm 



I. INTRODUCTION 

The detection of Dark Matter (DM) particles in ac- 
celerators or with direct and indirect searches would be 
of paramount importance for particle physics and cos- 
mology. A convincing identification of DM, however, 
will require a combined analysis of different experimen- 
tal strategies (see Refs. [1-4] for recent reviews). From 
LHC data alone, for instance, it will be difficult to recon- 
struct the relic density of the newly discovered particles 
(if any) , even if one assumes to know the underlying the- 
ory (e.g. Ref. [5]. See also [6] and references therein). 
However, combining LHC data e.g. with direct detection 
data significantly improves the reconstruction [7]. One 
could cite several other examples of complementarity of 
different detection techniques, such as the combination 
of different targets in direct detection experiments [8- 
11], or multi- wavelength analyses in indirect detection 
experiments [12]. 

Although the combination of different detection strate- 
gies provides an effective way of removing degeneracies in 
the parameter space of DM model, a combined analysis 
is made complicated by the fact that one needs to put 
together results that depend on different physical quan- 
tities, or on the same quantities but in a different way, 
a circumstance that makes the treatment of systematic 
errors (i.e. depending on the underlying model assump- 
tions) rather problematic. 

The statistical and systematic uncertainties affecting 
the reconstruction of DM particle from various experi- 
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ments have been discussed widely in literature. One re- 
cent example is the the uncertainty on the local DM den- 
sity, affecting both the rate of nuclear recoils in direct 
detection experiments, and the flux of upgoing muons 
induce by DM annihilations in the Sun. The statistical 
error on this quantity, based on dynamical constraints 
of various tracers in the Galaxy, have been discussed in 
Refs. [14-17] and more recently in [18, 19] (see also Ref. 
[20] ) . The systematic error due to the poor knowledge of 
the DM density profile in the Galaxy, and in particular 
to its possibly triaxiality has been discussed in Ref. [21]. 
Another recent example is the solar model dependence of 
the DM-induced neutrino signal [22]. 

Here we want to assess a minimal level of uncertainty in 
the normalization of the neutrino flux from the Sun (and, 
a fortiori, the Earth) due to the ignorance on crucial as- 
trophysical parameters, in particular concerning the DM 
phase space distribution. While the impact of these un- 
certainties for direct dark matter searches has been long 
acknowledged (see e.g. [23]) and it is still actively inves- 
tigated (see for example [24, 25]) the impact for indirect 
neutrino detection has been investigated less extensively. 
Occasionally, it has been noted that some effects (e.g. 
dumpiness and dark disk, see Sec. HI. 5 below) may en- 
hance the "baseline" flux, but little attention has been 
paid to the robustness of the baseline flux from an ordi- 
nary, smooth halo. Also, it is of interest to address the 
question of what would be the uncertainty on the pre- 
diction for the indirect neutrino signal, should a direct 
detection be used for normalization. Finally, since this 
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uncertainty turns out to be the dominant one, its value 
also sets the accuracy needed for several nuclear and par- 
ticle physics input parameters, derived either empirically 
or via theoretical computations. 

Here, we study this issue also in light of recent re- 
sults from numerical simulations. The structure of the 
article is the following: in Sec. II we review the depen- 
dence of the recoil and muon rates on key astrophysical 
parameters. In Sec. Ill we discuss the systematic un- 
certainties inevitably associated with these parameters, 
including time-dependent effects. Sections II and III are 
intended mostly as reviews of the essential concepts and 
quantities and to introduce our notation. In Sec. IV 
we discuss the impact of uncertainties on four variables 
discussed in Sec. Ill on theoretical predictions for the 
recoil and neutrino rates. We also illustrate the qualita- 
tive changes in the absolute and relative uncertainties as 
a function of the dark matter mass as well as direct de- 
tection experiment target material and energy-threshold. 
Finally, in Sec.V we discuss the results and present our 
conclusions. 



of gravitationally bound DM particles (see [28, 29] and 
refs. therein). 

We shall not enter the details of the form factor depen- 
dence to avoid introducing too many independent param- 
eters: for our purposes, we shall use the parameterization 
in terms of the nuclear mass from [1] 



Fl{E) = e^vi-rriA RaE/S) 



(3) 



with Ra = lfm[0.3 -t^ 0.91 ^mA,Gcv]- The over- 
all rate R of interest hero is obtained by integrating 
Eq. (1) over energy, starting from the threshold Ex, 
which is the minimum recoil-energy detectable, which is 
experiment-dependent but usually of the order of ~ 10 
keV. So, the input parameters to specify in this case are: 
{m^, E't, m^^, c7o,A, y}, where y denotes collectively the 
astrophysical parameters. 



Muon rate induced by neutrinos from DM 
annihilations in the Sun 



II. THEORETICAL SETUP 

The theory of DM experiments has been widely dis- 
cussed in literature. We summarize in this section the 
dependence of the predicted detection rates on some key 
parameters. 



As mentioned above, the detection of a high-energy 
neutrino flux from the center of the Sun or the Earth 
would provide a very convincing evidence for DM anni- 
hilations. The rate of events in neutrino telescopes de- 
pends on the number of particles that are captured, and 
on their fate inside the star or planet. 

The number of DM particles captured in the 
Sun/Earth obeys the following time evolution equation 



A. Recoil rate in direct detection experiments 

We start from direct detection experiments, for which, 
unless otherwise stated, we follow the notation of [26]. 
Apart for numerical constants of order one, the differ- 
ential number of recoil events per unit time and target 

mass writes 



dR 
dE 



dv^Qiv - v^i^) (1) 



where ij,^^a = MA/{m^ + Ma), Vmin = \J ^'^'^ 



h{v) = j dQt;V(v) 



JO 



dvh{v) = l, (2) 



where is the DM mass. Ma the target atom mass, 
(To, A the DM-nucleus cross-section. Fa the form factor, 
Pq the local density of DM, and /(v) is the 3D velocity 
distribution. Note that fi{v) denotes the velocity dis- 
tribution of DM parti(;lc!s at the Earth. If one has the 
distribution of velocity in the "halo rest frame" , in gen- 
eral one has to correct for: i) time-dependent velocity of 
the detector with respect to this halo, including Earth 
motion with respect to the Sun and Solar motion with 
respect to the halo; ii) Gravitational effects, including in 
general both the gravitational focusing of unboimd par- 
ticles (see e.g. [27] and refs therein) and the population 



N = C- CaN'^ 



(4) 



where C is the capture rate and Ca (see below) regu- 
lates DM annihilations, //both coefficients are constant, 
solving for N{t) one can derive 



^A{t) = ^tanh^ ^-^ J , ,eq 



(CCa)-'^'. (5) 



In the limit where steady state is reached within 
timescales much shorter than the lifetime of the Sun, 
^0 ~ 4.6 X 10^ yr, one has N = 0, N = y^C/CA, and 
the annihilation rate writes 



c 

2" 



(6) 



In this regime, the normalization of the signal only de- 
pends on C. More generally, the; present value r^(f0) 
depends also on Teq, i.e. on Ca- Ca can be written in 
terms of effective volumes V12, as 



Ca 



3 mpjTo 



3/2 



(7) 



with To , po respectively the central temperature and den- 
sity of the body under consideration [1]. Note that the 

above formula already has an uncertainty due to the only 
approximate assumption of thermalization, homogeneous 
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conditions of the core, etc., some of which are known to 
fail in some circumstances, see e.g. [30]. In the other 
useful limit Ten '> tp 



©I 



(8) 



cq 



This has the following implications for the large r, 
regime: i) the signal is quadratic in C; ii) there is an addi- 
tional dependence on Ca (i-C. linear in {(Tav){ ^"^'^° j"^^^) 
as well as a quadratic one on . Even neglecting particle 
physics unknowns, this means that additional astrophysi- 
cal imccrtainties enter. For the Sun, different models (see 
e.g. Rcf. [31, 32]) yield different predictions for these pa- 
rameters at the ^ 1% level. For the Earth, uncertainties 
are much larger. For example, according to [33], Earths 
temperature at the inner core boundary is estimated at 
5650 ± 600 K. We note incidentally that despite these un- 
certainties, it is possible to obtain useful constraints on 
DM particles with very large cross-sections, higher than 
those exchidcd by direct detection experiments, by re- 
quiring that they do not over-heat the Earth [34]. 

Focusing on C, we shall first remind that for most par- 
ticle physics models (including Kaluza-Klein and most 
neutralino DM models), the capture in the Sun is actu- 
ally dominated by spin-dependent interactions on hydro- 
gen [3]. This implies among others that we can consider 
the form factor « 1, which greatly simplifies the formu- 
lae. Although exceptions do exist, in the following we 
are interested in discussing the uncertainties induced by 
the astrophysical (rather than nuclear /particle physics) 
input, so we shall limit ourselves to this single contribu- 
tion in the capture. 

Prom rewriting Eqs. (2.8,2.13) of [35] in our notation 
we get 



^ — O^O.p X 



Jo Jo 



du 



1- 



where 



, / 4 my m-n , 



rriy — nir, 



.(9) 



(10) 



and u is the escape velocity from the unit shell volume 

considered, in turn depending on the distance from the 
center r; to a good approximation [13, 22], 



1355 km/s, Vs — 818km/s 



(11) 



where M.{r) is the mass within the radius r, in units 
of the solar mass. To derive the above equation for the 

capture rate, we replaced the Sun with a pure hydro- 
gen sphere and introduced the factor Cp < 1 to account 
for the fact that only a fraction of the Sun is made of 
hydrogen. What is crucial here is the approximation 
that the density of hydrogen traces the total one of the 
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FIG. 1: Timescale to reach equilibrium between capture and 
annihilation of DM particles in the Sun (Eq. (5)) as a function 
of the muon rate on a terrestrial detector for DM models in 
a pMSSM scenario, as implemented in DarkSUSY [36]. Solid 
circles correspond to models that fulfil the constraints on the 
cosmological relic abundance. The two horizontal lines corre- 
spond to the solar system age and 1/4 of the solar revolution 
time in the Galaxy. The vertical line shows the approxmatc 
sensitivity of IceCube plus DeepCore in 5 years of data taking 
[37]. 



Sun, which amounts to a small (^ 8%) overestimate of 
the capture [22]. Note that the consequences of our ap- 
proximations are further softened by the fact that we are 
only concerned with relative response of the observables 
to variations in the astrophysical inputs, so that overall 
scale mismatches are of no importance here. 

The main input parameters in this case are thus 
{m^, (To.pj y}- Further simplifications might be obtained 
by specifying the function /i, but we do not proceed fur- 
ther along this path, since we want to keep some gener- 
ality concerning /i, whose key parameters are described 
in the following section. 

Neutrino experiments such as IceCube are mostly sen- 
sitive to the fiux of up-going muons induced by high- 
energy neutrinos. The rate of /x's induced by a neutrino 
flux ^^{E^) is given by 



/ dE, 

Je}}^' Jo 



dyA{E^)P^{E,,y)^,{E^) (12) 



where y„ = 1 — E^^^/E^ and i?*'"' is the muon threshold 
energy of the experiment. A{Efj,) is the effective area 
of the detector and [Ei, , y) is the probability that a 
neutrino of energy E^ interacts with a nucleon producing 
a muon of energy Ef_t = {\ — y)E„ above the detector 
threshold energy, and can be estimated as 



R{E^,E';^^)a{E,,y) (13) 



where p is the density of the medium (typically wa- 
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ter/ice), R{E^, E^^') is the muon range in that medium, 
i.e. the distance traveled by muons before their energy 
drops below E^^^^ , and a{E^,y) = da'QQ{Ey,y)/Ay is the 
differential cross section for neutrino-nucleon charged- 
current scattering (for further details see e.g. Ref. [3] 
and references therein). In our case, the neutrino flux is 
simply given by 



47rD2 dE 



(14) 



where dN^/dE is the neutrino spectrum per annihilation 
and D is the distance of the Sun to the Earth (Astro- 
nomical Unit). 

We show in Fig. 1 the timescale to reach equilibrium 
between capture and annihilation as a function of the 
muon rate on a terrestrial detector for DM models in a 
pMSSM ( "Phenomenological Minimal Super- Symmetric 
Standard Model") scenario, as implemented in Dark- 
SUSY [36]. For reference, we also show with horizontal 
curves two characteristic timescales of the solar system: 
its age, and 1/4 of the revolution time around the Galac- 
tic center Tgai, i.e. the time to span a quadrant. The 
vertical line shows the approxmate sensitivity of IceCube 
plus DccpCorc in 5 years of data taking (e.g. Ref. [37]). 
It is clear that most interesting model in this scenario 
have Teq G (0.2; 1) Tgai, which makes "galactic" effects po- 
tentially important. Also note that for the Earth equi- 
librium is rarely reached in interesting models, see for 
example Fig. 18 in [28]. 



III. ASTROPHYSICAL UNCERTAINTIES 

In general, both DM velocity and density distributions 

crucially affect the the predictions for observables. Since 
in the literature most attention has been paid to the com- 
mon uncertainty coming from the local normalization 
value of the DM density, here we shall be mostly con- 
cerned with astrophysical quantities affecting differently 
direct and indirect signals. In particular, in the following 
we shall treat mostly velocity distribution uncertainties 
(subsections 111.1,2,3). In subsections 111.4,5, however, 
we shall also comment on the role of time-dependent ef- 
fects, which possibly introduce a "relative bias" between 
local density in direct experiments and long-term, aver- 
aged properties probed via capture. This is a subject 
rarely mentioned in the literature. 

Preliminarly to a more detailed discussion, it is worth 
recalling some facts about a crucial ingredient to make 
predictions for, or intepret, DM experiments: Vc{Rq), 
the velocity that a test-mass would have on a circular 
orbit in the Galactic Plane at the solar distance Rq from 
the GC ^. Note that, to a good approximation, the Milky 
Way mass distributions is axisymmetric: one can flnd a 



cylindrical coordinate system {R, 4>, z} centered on the 
GC where the potential is independent of with the 
z = plane dubbed "Galactic Plane", (see below the 
discussion of deviations from axial symmetry). For a 
particle having initial position at a distance R from the 
GC. with no radial and no vertical component of the ve- 
locity vector, the motion reduced to one in an attractive 
central potential. For any distribution of mass one can 
find a value of the azimuthal velocity Vc{R) which keeps 
the particle on a circular orbit at distance R, known as 
circular velocity. Formally, in analogy with the formula 
for spherically symmetric systems, one can always write 



GM _ v^{Rf 



R 



(15) 



where M is a parameter depending on R as well as the 
form of the potential. Actually, for axisymmetric poten- 
tials, the mass parameter M represents the mass enclosed 
within the radius R only for Mestel's disks, i.e. for disks 
whose surface density scales like S(i?) oc 1/R, while in 
general Vc will depend also on the mass outside R. For 
instance, in the case of an exponential disk, the above 
equation underestimates by up to 15% for radii larger 
than the scale radius of the disk [38] . Anyway, real stars 
have orbits varying in all coordinates. For stars hav- 
ing initial velocities in the plane z = 0, one has exactly 
(vz) = 0, while the orbit has a "rosetta" shape in the 
{R, (j)} plane, with R varying between a minimum and 
a maximum value. If the star moves off-plane, then an- 
other oscillation takes place above/below the plane (for 
more details, see e.g. Chap. 3 in [38]). So, in general to 
derive Vc requires identifying the circular motion of the 
"local standard of rest" (LSR), with respect to which ac- 
tual stars in our neighborood (including the Sun) possess 
relative motions. In [18], this was evaluated from a com- 
pilation of data to be Vc = 245 ± lOkm/s, but note that 
other observables may lead to a lower values, closer to 
220 ± 20km/s recommended by the lAU [39]. For the 
purposes of our estimate, we assume a central value of 
Vc = 235km/s, affected by a ~ 10% error. How reason- 
able this is can be also inferred from the collection of 
data shown e.g. in Ref. [40]. 



1. The halo typical velocity 

The halo velocity distribution cannot be directly ob- 
served, but must be inferred from a combination of obser- 
vations and theoretical considerations. In this inference, 
several sources of error enter. First of all, one assumes 



^ Actually, VcIRq) only matters with the Ansatz that the local 



phase-space density is dominated by a smooth halo component 
that is at rest in a Galactoccntric frame. Also implicit through- 
out this paper is that only isotropic velocity distributions are 
considered. Evidently, more realistic assumptions can only in- 
flate the error budgets with respect to those considered here. 
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that the rotation curve is constant over a significant ra- 
dius around the solar position, i.e. Vc — const. As shown 
by visual inspection of Fig. 1 in [40] , this approximation 
holds at not better than ~ 10% level. For a spherically 
symmetric system this implies that the gravitational po- 
tential $ is logarithmic, since 

v1{r) = r— => <l?(r) = v1 log(r) + const. (16) 

Until now, we assumed nowhere that the DM dominates 
the overall potential. If we did so, we would find the 
isothermal profile p oc r^^ as unique solution satisfy- 
ing the above (approximate) observational condition and 
symmetry condition, and the Maxwell-Boltzmann (MB) 

distribution /(v) a exp with vq = Vc as veloc- 

ity distribution. However, fit models suggest that DM 
only contributes a fraction of the v-^ at the solar system. 
For examples, the two models in table 2.3 of [38] yield a 
fraction of 38% and 64%. Note also that this fraction is 
changing with radius, too. So, in order to derive the link 
between, say, the average dispersion in DM velocity and 
the DM profile one needs further input, either from fits 
to the data or simulations. For a density of DM which 
goes as power-law p(r) oc r^^ in the neighborhood of 
th(^ Sun, the Eddington equation (see again [38]) yields 
a phase space distribution function for DM which has a 
dispersion 



(17) 



it is immediate to check that assuming that the only 
u— dependence enters as /(v) oc exp ^— ^j, one gets 



d3v^;V(v) 



2 



^^0 




(18) 



where the second equality is required for Eq. (17) to hold. 
Note that it has been often acritically assumed (perhaps 
following historical papers as [41] and [1]) that vq = Vc, 
but this strictly requires /3 = 2 at the solar radius. If 
one fits the DM profile with another popular model from 
numerical simulations, say NFW [42] with a scale radius 
~ 20kpc, there is a departure from this value of about 
20%. Alone, such uncertainty on the radial profile of 
DM at the solar radius implies an error of ^ 10% on 
the equality Vc = Vq. As described previously, a similar 
level of error appears inherent to the assumption of the 
constancy of Vc. Other errors, as the limitations of the 
spherical approximation (and thus isotropy of the space 
and velocity dependence of the DM profile) are more dif- 
ficult to assess, but likely present as well, with compara- 
ble or larger amplitude. We shall comment on that a bit 
further in the following. 

In summary, an estimate of the typical velocity param- 
eter vq is 



with a ^ 10% value due to the uncertainty on the mea- 
sured Vc (and thus correlated with its true value) and 
probably an uncertainty at least twice as big which is 
inherent to the theoretical assumptions. The latter er- 
rors are unfortunately of systematic nature and hardly 
reducible in the near future. 



2. Escape velocity 

From model- dependent considerations similar to the 
ones of the above section one can also derive that ex- 
pressions for the escape velocity as a function of Galactic 
parameters, see e.g. [41]. Equivalently, based on milder 
assumptions and on surveys of high velocity stars, one 
can infer constraints such as [43] 



498km/s < ■ymax(-R0) < 608km/s, 



(20) 



with a median value of 544km/s [43] which we adopt 
in the following as benchmark. Note that since theo- 
retically one expects that u^^x '^cj it is not surprising 
that a similar statistical error of ^ 10% is found for both 
Vmax and Vc- On the other hand, the local escape veloc- 
ity depends on global properties of the dark matter halo. 
Hence, to be conservative we shall treat the uncertainty 
on Vmax as independent from the uncertainty on vq as 
inferred from Vr- 



3. Shape of the velocity distribution 

In the above considerations, for simplicity we 
parametrized the uncertainty on the velocity distribu- 
tion via the typical velocity in the halo, or equivalently 
vq. But the event rates are not only sensitive to the 
second moment of the distribution, but also to its shape. 
There is no guarantee that the velocity distribution func- 
tion is truly MB, and actually several numerical simula- 
tions (see [25, 44]) suggest that noticeable deviations are 
present. To get a feeling for the sensitivity of the rel- 
evant observables to deviations from the MB shape, we 
parameterize the deviation from a MB. distribution as 



/(v,a) oc 



/mb(^^) (21) 



^^0 



(19) 



with a e [0.7; 1.3], as suggested e.g. from the right panel 
in Fig. 5 of [25]. Note that if the capture equilibration 
is long, the capture may probe the average distribution, 
while the direct detection probes the present one (see 
below) which is another effect linked to the possible shape 
deviation, not explored here. 

What we shall employ is the Eq. (21), with a MB cut 
at the iimax, correctly normalized to one. To that distri- 
bution one has then to apply a Galilean transformation 
to account for the motion of the Sun with respect to the 
DM halo, see below. 
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4- Motion of the Sun with respect to the DM halo and 
time- dependent effects 

Another effect on the recoil and capture rate is due to 
the fact that the Sun is in motion with respect to the 
DM halo "rest frame" with a velocity (since we are 
interested not in directional signals, only the modulus vq 
is going to matter in what follows.) Naively, one would 
expect 

Vq^Vc. (22) 

The situation is actually more complicated than that: 
the Sun is not rotating uniformly around the GC, i.e. it 
is not at rest with respect to the LSR frame. Then, one 
has a relative motion between sun and DM halo given by 

V0 = VLSR + Vloc . (23) 

From the rcc;cnt analysis in [45] one has the relative ve- 
locity determination Vioc = {U,V,W} = {11.1 ± 0.7 ± 
1, 12.2 ± 5 ± 2, 7.3 ± 0.4 ± 0.5}km/s respectively radi- 
ally inwards (U), in the direction of the Galactic rota- 
tion (y), and vertically upwards (W), with the former 
error statistical and the second one systematic. Taken 
at face value, this only amounts to an upward correction 
of about ~ 6% to the naive estimate Vq w w^, with an 
uncertainty which is even smaller, thus providing a sub- 
leading error. However, the situation is different when 
time-dependent effects are taken into account. In fact, 
as we have seen, one has to contemplate the possibility 
that the equilibration time for capture in the Sun is com- 
parable to (or longer than) a fraction of its orbital time 
^ 2 X 10^ yr. This situation is verified in a significant 
fraction of SUSY models leading to a large muon flux, 
as shown in Fig. 1 (see also the discussion in Ref. [22]). 
For the signal of the Earth, this is true for practically all 
interesting models. 

The fact that the capture signal depends on the long- 
term past history of the Sun (and Earth) immediately 
changes the importance of the solar motion effect men- 
tioned previously: differently from the first term at the 
RHS of Eq. (23), the second term is time-dependent due 
to the variation of velocity of the Sun along its orbit in 
the Galaxy, with typical timescalc of ^ 2 x 10^ yr. 

To be more quantitative, one would need to know the 
actual orbit of the Sun in the Galaxy and the equilibra- 
tion timescalc for the case realized in nature. Both have 
significant sources of error: i) the error in the "initial" 
(actually final, since the evolution is performed back- 
wards!) conditions, in particular the velocity vector Viod 
ii) the only approximate knowledge of the Galactic po- 
tential: iii) the incomplete knowledge of particle physics 
and astrophysical parameters determining Teq. As a re- 
sult, rather than "a refinement" to the naive estimate 
we can conservatively consider this kind of effects as an 
additional source of error. For illustration, we take the 
results of the numerical integration performed in [46] for 
the model 2 in [47] (note that compared with those re- 
sults the local velocity component has been significantly 



revised in the analysis reported in [45]). As a result, the 
distance from the GC varies between -2% and +7% from 
the present value, while the velocity results to anticorre- 
late almost exactly, with an equal (in modulus) and oppo- 
site sign variation (the motion is indeed quasi- keplerian). 
As a consequence, the two effects almost c;ancc!l out for 
the capture rate in the Sun, so this time-dependence is 
not the main source of error within the simple model con- 
sidered. 

The situation is actually different if we move beyond 
a smooth and isotropic halo approximation: dropping 
the isotropy means that the local density might differ 
from the averaged one (to which the Solar capture is 
sensitive for long equilibration times). This difference 
has been recently evaluated in [21] based on numerical 
simulations: in the ones including approximatively the 
effects of baryons, this error is about 20% to 30%. In 
those without baryons, depending on the ignorance of 
our orbit within the triaxial halo of DM, the error is 
as large as a factor ^2. To account for triaxiality, a 
simple model with a dependence of the density along the 
Galactic Plane as p{(f)) = (/o)gai[l + £ sin(27r(/))] can be 
taken as representative of results with baryons, for e = 
0.3. Of course, since we do not know the phase at which 
we are at current time, one can define an effective density 

(p)eff(0) = ^r p(<A + ^)d^ (24) 

which is probed by the capture rate and is a function of 
(f). We show in figure 2 the effective density as a func- 
tion of position in a circle of size Rq around the GC, 
for different choices of Teq in units of Tgai. We also show 
in the same figure (right panel) the maximum relative 
difference between the effective and instantaneous den- 
sity, as a function of the time to reach equilibrium Teq, 
in units of a galactic year Tgai. This is obtained at each 
point by scanning over the unknown (/). In both cases, we 
fix £ = 0.3. Two comments are in order: the maximal er- 
ror obtained when Teq/rgai is very large reflects the simple 
fact that the capture probes the real average value (p)gai, 
while the recoil the instantaneous value which can differ 
as much as e = 0.3 from the average (black, solid curve 
in the left panel of Fig. 2). However, the situation is even 
more interesting when Teq/rgai ^ 0.6 (which Fig. 1 hints 
to be a physically interesting range): in this case the rela- 
tive error can range from almost negligible to above 40%, 
within only a factor 3 variation of Teq. This suggests that, 
even if we knew that the ellipticity of the Milky Way halo 
is at the 30% level, and even by assuming that an even- 
tual neutrino detection hints at Toq/Tgai ^ 1, the lack of 
a detailed knowledge of the particle physics and astro- 
physics limit the estimate of this error conservatively at 
the 40% level. Needless to say, should the baryon disk 
effect be less prominent than what suggested by simula- 
tions, this effect could be a factor of several larger. The 
actual impact of baryons on the DM structure is in fact 
still subject to debate. 
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FIG. 2: Left Panel: Effective density as a function of the position in a circle of radius 7?© around the GC, for different choices 
of Teq . Right Panel: Maximum relative difference between the effective and instantaneous density, as a function of the time to 
reach equilibrium Teq, in units of a galactic year Tgai- 



5. Additional halo (suh) structures 

Concerning the non-smoothness of the halo, an at- 
tempt to address its impact on detection rates was per- 
formed in Ref. [48], where the effects of the enhanced 
density while crossing a sub-halo were considered. For 
our purposes of assessing reasonable errors, rather than 
exploring what can happen in principle if such an en- 
counter takes place, it is more interesting to assess the 
probability that such an encounter takes place at all. If 
we focus on subhalos of mass = 10~^Mq, we can 
work out the characteristic size, number density there- 
fore average separation for these objects. The Sun en- 
counters a subhalo of mass M^h and size with a fre- 
quency Fsh = 47rutot'T-sh?'sh- Assuming that a fraction / 
of the DM is in the form of clumps, and inserting typical 
values, we find 



220km s-1 Msh VO-Olpc; ^ ' 



To perform this estimate righ we have followed the usual 
strategy of extrapolating the mass function of subhalos to 
very small masses, and assumed that the spatial distribu- 
tion of subhalos traces the smooth DM distribution. The 
second assumption is most likely a very crude approx- 
imation, as the most recent numerical simulation have 
shown that the distribution of subhalos is anti-biased, 
such that if p{r) is the total (smooth plus subhalos) DM 
density profile, Ushir) oc rp{r). The fraction / of DM 
in subhalos is therefore a function of r and it tends to 
deplete the number of clumps at the solar radius. The 



results of Ref. [49] point towards a more realistic value^ of 
/ ~ 10~^. This estimate suggests that the timescale for 
an encounter is comparable to (a relevant fraction of) the 
orbital time of the Sun in the Galaxy. On the other hand, 
the crossing time tsh for such an encounter is very short, 
of the order of a century. Hence, already on this basis we 
can conclude that significant alterations of the annihila- 
tion yield, while possible in principle, are very unlikely. 
In the above considerations, we implicitly assumed that 
the time during which a significant enhancement of the 
signal takes place is comparable with the crossing time. 
In Appendix A we show that this is a good assumption if 
large boosts are required, since a decay time much longer 
than the crossing time is only possible at the expense of 
reducing the enhancement in the signal by the same fac- 
tor. 

We also note that the above estimate is not very de- 
pendent on the exact mass of the sub-halo considered: 

1 /3 

since one has roughly rgh ~ , and tgh oc rgh, it fol- 

— 1/3 

lows the weak dependence Fgh ~ , while Fgh ^sh 

which roughly quantifies the probability of "living dur- 
ing a crossing" stays virtually constant and very small, 
of the order of 10~^ for the parameters quoted above. 

It is worth mentioning that it is not even sufficient that 
an encounter takes place in order for a significantly en- 
hanced capture to happen. In fact, as a first approxima- 
tion we can think of each subhalo as orbiting the galaxy 



^ It is perhaps curious to note that if an unrealistically large value 
of / is used, the actual annihilation signal could be lowered with 
respect to its naive equiUbrium value in the limit of smooth 
halo [48]. 
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with a typical velocity w sampled from the /(v) distri- 
bution. Neglecting the Wmax cutoff and assuming a MB 
distribution, the ID distribution at the Sun is readily 
written as 



2w 



/IT Vo Vq 



■ sinh 



2wvq 






exp ^ 


V '^0 J 





(26) 

The particles in each sub-halo have a dispersion a around 
w. However, as long as cr ^ wq (as expected for the small- 
est, hence more abundant, sub-halos) this only brings in 
a small correction. Roughly speaking, as long as long 
w < Umax the DM particles will be captured, when in- 
stead in an encounter one has w > Umax only a very 
small fraction of the DM population could be captured 
(see Eq. (10)). For fiducial parameters in Eq. (26) and 
taking into account that for the Sun Umax — 200km/s 
for a 100 GeV particle, one finds that only ^ 15% of 
the encounters will lead to a(n efficient) capture. For the 
Earth, this fraction is negligible. 

The scenario could be further enriched if one consider 
the possible existence of the "dark disk" [50, 51] due to 
the dragging effect of the baryonic disk on the DM halo. 
The crucial aspect is that the (typically subleading by 
number) population of DM disk particles has a relatively 
low lag velocity with respect to the stars, uiag ^ — 
150 km/s, and a comparably low velocity dispersion a^. 
As a consequence, the Sun and the Earth effectively see 
a colder DM gas, a fact that eases captures. For the 
Sun, the scaling relation for the capture reported in [51] 
suggests that, depending on parameters, the annihilation 
signal can obtain a correction which ranges from marginal 
(a few percent) up to more than one order of magnitude. 
For the Earth the dependence is by far more dramatic, 
usually amounting to orders of magnitude. This confirms 
once again that the normalization of the signal from the 
Earth is extremely dependent on the details of the (very) 
cold tail of DM phase space distribution and, while still 
potentially useful for a serendipitous discovery of DM, 
can tell us very little on the particle physics. 



IV. HALO UNCERTAINTIES ON DM 
DETECTION RATES 

In most of this section, we estimate the variations in 
the direct experiment recoil rate R as well as capture 
rate C when velocity distribution parameters vary within 
their uncertainties, as discussed in Sec. 111.1,2,3. At 
the end of this section, we come back to the additional 
uncertainties induced by time-dependent effects on the 
effective density probed, in particular those mentioned 
in Sec. III.4. 

The two functions R and C depend in general from 
non-astrophysical parameters x as well as some as- 
trophysical parameters y. Following the discussion 
in Sec. II, we select for the former parameters x = 
{ui^jEtjA} denoting respectively: i) the DM mass, 
which enters both R and C, whose fiducial value is put 
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FIG. 3: Impact of varying simultaneously the four astrophys- 
ical parameters y(see text for range of variation), in the R — C 
plane, for different values of the DM mass (top) different val- 
ues of the energy threshold (center) and different values of 
target mass A (bottom). In each panel, the unmentioned pa- 
rameters are held fixed at their fiducial values, = 100 GeV, 
Eth = lOkeV, and A = 73. 
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to 100 GeV. ii) two parameters for the direct detection 
experiments, namely the threshold energy i?th for de- 
tectable recoils and the atomic number A of the tar- 
get material. Both are model-dependent, but the for- 
mer is usually in the 10 keV range. We shall assume 
Eth = lOkeV as fiducial value, but we shall explore the 
consequences for R of varying it by a factor 2, upwards 
or downwards. The target material atomic number A en- 
ters our formulae for C via the form-factor dependence. 
We shall take the fiducial value A = 73 for an "inter- 
mediate mass" target, such as Germanium; A = 131 
for a heavy target material as Xenon; A = 28 for a 
relatively light target, such as Silicon (these identifica- 
tions should not be taken too seriously but to investigate 
the dependence on the target mass, since the form-factor 
functions are only schematic). Note that the DM-direct 
detector target cross section (typically spin-independent) 
cro,A and DM-solar matter target cross section (typically 
spin-dependent on hydrogen) (Jq.p enter just as overall 
normalizations in R and C; thus, exactly like the value 
of the local density of DM pq, they do not afi'ect our 
considerations. 

Following the discussion in Sec. Ill, we identify four 
main (time-independent) astrophysical variables y = 
{vo, ^'max, vq, a}, denoting respectively the halo disper- 
sion velocity, the escape velocity, the local velocity of 
the sun in the DM halo, and the shape-parameter of the 
velocity distribution. They are varied as follows: vq at 
fiducial value of 235 km/s as well as varied by ±10%, 
Vq at fiducial value of 235 km/s and with ±10% varia- 
tions, a at the values 0.7, 1, 1.3 and Wmax is set at 498, 
544, and 608 km/s, with the central values of the latter 
two considered as fiducial ones. Note that in particular 
the assumed 10% uncertainty on vq is a very conserva- 
tive estimate of the error on this quantity, with a more 
realistic value even a factor of few higher, as discussed in 
Sec. III.l. This is however consistently with our goal of 
assessing a minimal uncertainty on the observables. 

For a first glimpse at the overall dependence of the 
observables on parameters, we show in Fig. 3 a series 
of "scatter plots" in the i? vs. C plane scanning over 
the min-med-max values of the y parameters, which are 
varied simultaneously. These plots arc presented for the 
fiducial values of all but one of the x- variables: In par- 
ticular, the X variable varied in the top panel is the DM 
mass, in the central panel it is the energy threshold for 
direct detection, and in the bottom panel the target mass 
A for the direct detection. The observables are normal- 
ized to the values assumed for fiducial parameters. It 
is obvious that, in all cases, the uncertainty on the cap- 
ture can have excursions of ^ 40%, while the range span 
by R is typically within ±20% of its fidiicial value, al- 
though it can increase for light targets, high thresholds, 
or heavy particles. What is more important, there are 
(anti) correlations between variations in R and C, whose 
magnitude and even sign depend on the variable con- 
sidered. This is illustrated by diagonal "trends" in the 
dispersion plots, whose slope can change both quantita- 



tively and in sign. It is not unusiial that the relative 
uncertainty of the ratio C/R, say, is larger than the er- 
ror on each quantity C, R due to this effect. Just for an 
illustration of this effect, we computed the correlation co- 
efficient corr of the R 's and C"s for the scanned points. 
This is defined as follows 



corr{R, C) 



T.,(^i - m(ri - in) 



(27) 



where the index i denotes a particular realization of the 
arguments, and the sum runs over the the scanned pa- 
rameters sets. 

The dateset cross-correlation when varying m;^=20, 
100, 500 GcV varies from -0.12 to -0.15 and to +0.55; 
when varying i?th=5, 10, 20 keV changes from 0.003 
to —0.15 and to —0.42; when varying A from 28 to 73 
and 131 rises respectively from —0.50 to —0.15 and to 
±0.37. Of course, this has no truly statistical meaning 
(we assumed "flat" priors for the variables in their range 
and did not account for the existing correlations among 
them), but it clearly illustrates the trend of qualitative 
change of the correlation when underlying parameters are 
changed. 

For a more compact and quantitative exploration of 

the parameter dependence, it turns useful to define the 
logarithmic derivatives: 



7i(x, yo) 



51nC(x,y) 



(9 In Hi 



y=yo 



_ aini?(x,y) 

^i(x,yo) = 



dlnyi 



(28) 
(29) 



The functions 7^ and Qi are sufficient (to leading order) 
in describing the sensitivities to the quantities of inter- 
est. Note also that the change in C/R as a consequence 
of a change Si = Syi/yo can be written, in the linear 
approximation, as (no sum over repeated indices) 



dyi R 



li - Qi 



S,. 



(30) 



The values of the 7j and Qi coefficients for fiducial values 
of y and different values of the x are reported in Tab. I 
and Tab. II. 

From these tables it is clear that the dependence on the 
parameters is actually different in the two cases. Take the 
second row of Tab. I: it means that an uncertainty in the 
velocity of the Sun in the halo is responsible for a fac- 
tor I — 1.64 — 0.15 1 ~ 1.8 times larger uncertainty in the 
relative normalization of thc^ two rates. From that table, 
one immediately notices the general trend that, when the 
recoil rate only probes the high-energy tail of the distri- 
bution (light DM particles, high thresholds, light targets) 
the dependencies on virtually all parameters are opposite 
for R and C: the capture in fact always prefers the low- 
energy tail of the distribution, where particles are slower 
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-fi (20) 


g^ (20) 


-^^ (100) 


g^ (100) 


"f^ (500) 


g^ (500) 


Vo 


-0.80 


1.04 


-0.98 


-0.04 


-0.90 


-0.22 


Vq 


-0.94 


1.00 


-1.64 


0.15 


-1.92 


-0.05 


i'max 


-0.12 


3.21 


-0.13 


0.77 


-0.13 


0.50 


a 


-0.15 


0.22 


-0.16 


0.004 


-0.12 


-0.04 



TABLE I: The response functions for capture and direct re- 
coil rate with respect to different halo parameters for three 
values of the mass nix = 20, 100 , 500 GeV, and for the fidu- 
cial values of the parameters a = 1, ^ = 73, Eth = lOkeV, 
Vc = 235km/s, vq — 235km/s, Umax = 544km/s. Note that 
the product a p enters only as a scaling, and similarly for both 
observables (i.e. response is 1). 

and thus easier to capture. For higher target masses, a 
competing effect arises due to the form-factor suppres- 
sion, so that mixed trends arise. On the other hand, in 
the limit Eth — ^ 0, the dependence on the halo parame- 
ters goes to zero since the whole phase space is actually 
probed via direct direction. Note that the effect of the 
shape parameter a is rather modest, and the relative un- 
certainties due to the non-MB nature of the distribution 
can be estimated to be < 10%. A posteriori, this justifies 
a simple parametric approach to explore this variable. 



i 


g^{A = 28) 
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0.14 


-0.19 


-0.10 


0.11 


Vq 


0.32 


0.01 


0.07 


0.31 


Umax 


1.12 


0.52 


0.69 


0.95 


a 


0.05 


-0.02 


-0.01 


0.04 



TABLE IL As in Tab. I, for the direct detection rate, keeping 
rux = lOOGeV and varying the target mass and the threshold. 

On the top of the different sensitivity to the velocity- 
distribution variables, as already argued in Sec. IIL4, 
another source of uncertainty is the fact that the two sig- 
nals sample different time averages of the physical quan- 
tities, in particular the DM density. Accounting for a 
40% maximal error on the ratio of observables estimated 
on the basis of simulation results (see Fig. 2), we can 
conservatively conclude that there is likely a relative un- 
certainty on C/B. which reaches a factor of ~ 2, just due 
to "halo astrophysics" . 



V. CONCLUSIONS 

A conservative estimate of the error on the capture 
rate suggests that the normalization of the neutrino sig- 
nal from DM annihilation in the Sun/Earth brings not 
only the uncertainty given by the local DM density Pq, 
but often more importantly uncertainties connected with 
the velocity distribution function, as well as our motion 
in the halo. We showed that the error on the capture C 
can easily reach « 40%, which translates into a compara- 
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ble or larger error on the annihilation signal. Even more 
important, the relative uncertainty on the normalization 
of the annihilation to recoil signal is typically larger, eas- 
ily up to a factor two in presence of typical equilibration 
times of ~ 10'^ - 10^ years . For cases where equilibra- 
tion times are billions of years or larger (as typical for 
the Earth) the uncertainties are likely at least few times 
larger. The same happens if the effects of baryons in 
reducing Galactic halo triaxiality are smaller than what 
estimated at present from numerical simulations. 

We also studied the sensitivity of the observables to 
different input parameters, which was shown to be often 
opposite for the direct and indirect signals, due to the 
preferential probe of the high tail part of the velocity 
distribution for the recoil rate. Interestingly, at least for 
some variables this anticorrelation can be reversed de- 
pending on masses of DM and instrumental parameters. 

It is worth noting that these "halo" uncertainties ap- 
pear larger than uncertainties coming from solar compo- 
sition and/or nuclear /particle physics (for a recent esti- 
mate sec [22]) and thus provide in most cases the dom- 
inant limiting factor in the extracting particle physics 
information from a future detection of a DM neutrino 
signal, even when normalizing the rate to some direct 
detection event rates. This is true in particular if no 
or little spectral information is available. Enhancements 
in the annihilation rate are a generic expectation of the 
presence of a "dark disk" created by interaction of the 
halo particles with the baryonic disk, but its effect can 
range from few percent to orders of magnitude. The role 
of substructures appears instead marginal at best, bar- 
ring for highly unlikely circumstances. In principle, the 
velocity distribution of WIMP DM can be reconstructed 
from direct dark matter detection data in the range of 
velocities probed (provided the DM mass is sufficiently 
constrained), but O(IOO) events are needed to start im- 
proving the knowledge of fi(v) to better than what in- 
ferred indirectly via considerations as those reported in 
Sec. Ill (see e.g. [53]). 

The considerations developed in this article have sev- 
eral implications, whose exploration is left for future 
works. One obvious consequence is that great care must 
be paid when comparing the constraints on the spin- 
dependent scattering cross section inferred from direct 
detection with those inferred by neutrino experiments 
(see e.g. [54- 56]), especially when light DM particles are 
considered, as well as combinations of experiments with 
different target masses. This is notably the case of the 
comparison of the DAMA and CoGcnt results with the 
exclusion plots obtained by other experiments (see e.g. 
Ref. [57] and references therein). 

Another improvement over the first estimate provided 



^ Going beyond the WIMP paradigm and introducing quite large 
self-interactions between DM particles can cause a similar mis- 
match between the two types of observables [52]. 
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in this article would consist in providing a more statisti- 
cally sound assessment of the error range in the parame- 
ters governing the velocity structure of the halo (and the 
motion of the Sun through it), along the lines of what 
done in [18] for the local halo density. Here we implicitly 
assumed flat priors in the parameter space and uncorre- 
lated variations of the different parameters. Needless to 
say, such refinements would become compelling for ex- 
tracting particle physics parameters if the present gener- 
ation of WIMP DM detectors were to show any evidence 
for these particles. 
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Appendix A: Decay time vs. Crossing time 

In Sec. III. 5 we have assumed the crossing time of 
a sub-halo as the relevant timescale for enhancement in 
the annihilation signal. In principle, however, the rele- 
vant quantity is not the time spent by the Sun sitting in 
a dark-matter subhalo (crossing time), but the time for 
which the annihilation rate in the Sun is different from 
the case of a completely smooth halo (enhanced signal 
decay time). We show here that the latter can only ex- 
ceed the former at the expense of suppressing the signal 
enhancement. 

Let us assume that the Sun/Earth encounters a sub- 
structure with density K Pq, with K ^ 1, during a cross- 
ing of time tsh starting att = 0. The only case of interest 
is when a significant capture happens in the passage; in 
this case, the value reached by A'' can be significantly 
larger than its equilibrium (or long-term) value in the 
smooth halo, which we can neglect (i.e. N{0) « 0). 
Hence, in terms of the capture rate C in the smooth halo, 
the approximate solution writes (req = (C Ca)~^^^) 



N{t) 




^ tanh ( 



(Al) 



^eq 



As a consequence, the enhancement in the annihilation 

signal at the end of the crossing over the "naive" , long- 
term average signal in the smooth halo is given by the 
following "boost" function 



S(4h) = 



rA(^sh) 

ro(i0) 



tanh 



K- 



tanh' 



(A2) 



where the equilibration time is the one defined with re- 
spect to the smooth halo. 

For t > tsh, the enhanced capture rate ends, and the 
enhanced signal starts to decline. Its time evolution is 



now dictated by 



hence, for t > t^h, 



N 



-CaN^ 



(A3) 



7V(U) - N{t) = N{tsh)N{t) CA{t - U) • (A4) 

The signal (proportional to N^) will drop to a factor 
1/F of the one at tgh after a time delay Apt given by 



Apt = 



y/F-1 



(x/F-l)Te, 



CaA^(U) V^tanhfii^ 



(A5) 



Ideally, to maximize both the boost and the duration of 
the enhanced signal one would require that both B{tsh) 
and A pt are as large as possible (note that for any prac- 
tical circumstance (v^ — 1) ~ C'(l))- Clearly, these 
two requests are in tension with each other. There are 
two limiting situations (for simplicity, we assume in the 
following that the smooth halo signal has reached equi- 
librium): 

For equilibrium to be reached during the crossing, 
which guarantees the "full boost" , we require 



> 



' eq 



hence 



B{tsh)^K, Apt: 



(v^-l)r. 



eq 



K 



tsh ■ 



(A6) 



(A7) 



So, the optimal case for enhanced signal requires a "decay 
time" after crossing which is comparable or more rapid 
than the crossing time. 

If instead equilibrium is not reached, assuming the ar- 
gument of the 'tanh' function to be small one obtains 



Bitsh) = Ke, e«Jf( — ) <1, (A8) 

. To, 



'eq 



Apt « tsh ■ 



More in general, one has 



(A9) 



B{tsh)^ « (x/F- i)k'-^ < (^/F- 1)K (AlO) 



tsh 



y 



where y = y/K tsh/Teq- The inequality is saturated only 
for small y, while the expression at the LHS is suprressed 
as 1/y for large y. So, we conclude that a decay time 
much longer than the crossing time is only possible at 
the expense of reducing the enhancement in the signal 
by the same factor. 
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